Synchronization in Complex Networks: a Reply on a recent Comment 



S. Boccaletti 1 , M. Chavez 2 , A. Amann 3 and D.-U. Hwang 1 
1 CNR-Istituto dei Sistemi Complessi, 
Largo E. Fermi, 6, Florence, Italy 
2 Laboratoire de Neurosciences Cogmtives et Imagerie Cerebrale (LENA) CNRS UPR-640 
Hopital de la Salpetriere. ^7 Bd. de I'Hopital, 
75651 Pans CEDEX 13, France 
3 Tyndall National Institute, Cork, Ireland 

We clarify a number of points raised in [Matias, arXiv:cond- mat/0507471 v2 (2005)]. 
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In a previous comment on cond-mat 0, two of our 
recent publications 0, were compared with researches 
simultaneously carried out by the Potsdam group and 
published in 0, ■ It was argued that the conclusions 
of both sets of papers are closely related. To obviate 
confusion, we will in the following point out a number of 
facts which will help to better understand the differences 
in the results and conclusions of the involved works. 

1 . The root for the improvement observed in 0, is 
the normalization of the diagonal terms of the con- 
nectivity matrix. Such a normalization is limited 
in 0, 01 to the optimal case (3 = 1 , whereas in our 
papers it is realized for all values o/a | and 6 
Furthermore, the optimal condition [3 = 1 in 
corresponds to a = in and to 9 = in 
Therefore, any further synchronizability enhance- 
ment observed, e.g., in Figs. 2,3 of and in Figs. 
1,2 of cannot be related to any of the mecha- 
nisms reported in . It is therefore not appro- 
priate to compare the enhancement of 00 with re- 
spect to a non weighted configuration to the further 
improvements that our original approaches provide. 
For instance, the very nice analytical study of 00 
(duly quoted by us in 0) is restricted to a series of 
explicit hypotheses that require a sufficiently "ran- 
dom" configuration to be satisfied. Our studies 
manifestly show (see e.g. Fig. 1 and 2c of 0) that 
while for Erdos-Rcnyi graphs the maximal synchro- 
nizability is attained in the condition of (that 
is for 8 — 0), the interplay between hierarchy and 
weighting is able to further improve this value for 
other network's configurations with the same prop- 
erty of constant total strength of input connections. 
This extra enhancement cannot be extracted by the 
arguments raised in 0, . 

2. A weighted wiring does not necessarily imply asym- 
metric connections between elements. For instance, 
the weighting procedure of 0, 0] (based on powers 
of the node degrees fcf ) yields always (i.e. for all f3) 
a symmetric connection between connectivity peers 
(as e.g. it can frequently occur in assortative net- 
works). At variance, in [3j we are able to modulate 



and direct the interaction between any pair of nodes 
by the parameter 9 G ( — 1,1) (from a unidirectional 
to a bidirectional symmetrical one), accordingly to 
the hierarchical (age) ordering of the nodes. The 
conclusion that a weighted configuration based on a 
hierarchical structure can enhance synchronization 
is, therefore, original. 

3. Our studies 0] report a detailed analysis of the 
propensity for synchronization in parameter spaces. 
Namely, in (0) , we detail in Fig. 2 (Figs. 1,2) 
the behavior of Aat/A 2 (A^/Aj) in the parameter 
space [a, B] (#), showing that the enhancement due 
to a weighting procedure based on the link loads 
(the age ordering) is a generic feature in parameter 
space. In contrast to what was argued in 0, we 
in fact carefully considered the dependence of our 
results on the network size N in both our papers. In 
particular, we show via Gershgorin's circle theorem 
that the normalization procedure yields an upper 
bound for Aat for any network size and topology. 
Notice that, in 00, when the weight is selected 
to be fcf , A at is in general size dependent for j3 < 0, 
and one has limjv^oo A at = oo. 

4. The Author of is incorrect when stating that the 
main conclusion of is that the enhancement in 
synchronization is evident when the coupling direc- 
tion is from the younger to the older nodes. Our 
Letter explicitly claims that one has to combine 
this property (relaxed also in for (3 > 0) 
with a structure of interconnected hubs induced by 
the hierarchical (age) ordering. The second part 
of our Letter (as well as Fig. 2c of 0), is en- 
tirely dedicated to study a series of ad-hoc mod- 
ified networks with the aim of pointing out how 
such an interplay works in scale- free networks. Fur- 
thermore, it is straightforward to show that for 
any network size and topology one can construct 
a coupling scheme based solely on age ordering and 
able to give = 1 (the minimal as possible 
value). This can be realized by constructing a uni- 
directional tree structure, by means of which one 



2 



randomly selected node forces a cascade of nodes 
spanning the whole networks. Using properly nor- 
malized weights for each uni-directional connection 
(i.e. normalizing the connection strength to the 
in-degree of the slave node), and imposing a zero- 
row sum configuration induces a triangular matrix 

whose diagonal elements are (0,1,1,1,1,1,1, ,1), 

giving the following set of eigenvalues Ai = 0, 
A, = l,i = 2,...7V Q. 

5. It is well known that a crucial difference in topo- 
logical information content exists between the node 
degree and the load of the link or edge betweennes, 
the former retaining only information on the lo- 
cal structure of the network, the latter providing a 
measure of centrality of the node Q in the global 
topological structure. It is also well known (e.g. in 
the study of network's community structures) that 
nodes with high degrees do not necessarily have 
links with high loads, and links with high loads 
do not necessarily connect two nodes with high de- 
grees. An important and original result of our stud- 
ies is therefore that (under the assumption of 
input normalization) the knowledge of global topo- 
logical properties of the network can be used to 
improve the propensity of synchronization over the 
limit obtained in jj, |£j , where only local properties 
were used. This improvement is not closely related 
to the ideas or conclusions of 0, 0| . 

In summary we have shown that the results and con- 
clusions in Refs. 0, El 0, are independent and com- 
plementary. However, though from a completely differ- 
ent point of view than that suggested by [l| , a relation- 
ship between these works indeed exists. Taken together, 



the approaches of 0. 0. [J. Q have undoubtedly and in- 
dependently established weighted networks as a promis- 
ing framework for the further study of synchronized be- 
havior. While the arguments provided in the present 
Manuscript demonstrate that it is erroneous trying to 
deduce the conclusions of one approach from the other, 
we believe that much effort will be stimulated in the 
forthcoming years to investigate other ad-hoc weighted 
configurations, for assessing more deeply the underlying 
mechanisms for the formation of collective dynamics in 
all those circumstances where the weighting features can 
be directly extracted from available data. 

This work was supported by MIUR-FIRB project n. 
RBNE01CW3M-001. 



[1] M. Matias, |arXiv:cond- mat /050 747i"V 2 (2005). 

[2] M. Chavez, D.-U. Hwang, A. Amann, H.G.E. Hentschel 

and S. Boccaletti, Phys. Rev. Lett. 94, 218701 (2005). 
[3] D.-U. Hwang, M. Chavez, A. Amann and S. 

Boccaletti, Phys. Rev . Lett. 94, 138701 (2005); 
|arXiv:cond-mat /0412728| 
[4] A.E. Motter, C.S. ZlTou and J. Kurths, Europhys. Lett. 

69, 334 (2005); arXiv:cond-mat/0406207 
[5] A.E. Motter, c's. Zhou and J. Kurtns7"P hys. Rev. E 71, 

016116 (2005); |arXiv:cond-m at/0502309 
[6] The full description of this latter weighting procedure will 

be reported by us in a forthcoming paper. 
[7] L. C. Freeman, Social Networks 1 215 (1979); M. E. J. 

Newman, Phys. Rev. E64, 016132 (2001); U. Brandes, 

Journ. of Mathematical Sociology 25, 163 (2001); K. I. 

Goh, B. Kahng, and D. Kim, Phys. Rev. Lett. 87, 278701 

(2001). 



